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SUMMARY  &  CONCLUSIONS 

In  order  to  understand  the  problems  intractable  problems 
nature  presents  us,  we  are  forced  to  make  simplifications  till 
we  arrive  at  simple  analytical  models  that  we  are  capable 
of  understanding.  Such  canonical  models  are  rare,  but  useful 
as  tools  for  constructing  more  realistic  models  that  we  can 
use  analyze  nature.  The  class  of  analytical  models  for  clutter 
analysis  limited  to  those  that  consist  of  various  amplitude 
models  with  the  phase  noise  assumed  to  have  a  probability 
density  function  that  is  uniformly  distributed.  These  analytical 
models  can  be  extended  by  relaxing  the  assumption  of  uniform 
phase  noise  to  phase  noise  with  non-uniform  distributions.  It 
is  shown  how  to  determine  the  probability  density  function 
for  these  non-uniform  distributions  in  general  and  then  the 
method  is  illustrated  with  such  distributions  as  Gaussian, 
Laplacian,  and  Chi-squared.  This  enables  one  to  determine 
the  probability  density  functions  of  the  individual  components 
of  clutter  model  such  as  xt  =  it  cos(</>c). 

Using  the  rule  for  determining  the  products  of  distributions, 
we  show  that  the  functional  form  for  xt  is  reduced  to  eval¬ 
uating  integrals  that  reduce  to  elliptical  functions.  Once  the 
functional  form  has  been  determined,  it  is  easy  to  determine 
the  moments  of  the  PDF  and  hence  completely  characterize 
its  statistics. 

1.  INTRODUCTION 

Analytical  clutter  models  have  proven  useful  in  both  the 
design  and  analysis  of  moving  target  indicators  (MTI)  and 
pulse  Doppler  (PD)  radars.  Usage  of  them  to  represent  ran¬ 
dom  clutter  enables  one  to  predict  some  aspects  of  system 
performance  analytically,  which  is  always  useful  as  means  of 
testing  simulations  and  to  provide  sanity  checks  for  ones  data 
[1],  [2],  It  is  not  generally  realized  that  most  clutter  models 
can  be  cast  in  the  form 

ct  =  vt  cos(wcf  +  </>c)  (1) 

where  "hat"  denotes  random  variable,  Ct  =  distributed  clutter 
process  that  is  assumed  to  be  strict  sense  stationary,  i>t  A 
clutter  envelope  process,  loc  =  clutter  frequency,  (j>c  =  phase 
noise  process  (which  we  assume  is  not  uniform).  (I  will  refer 


to  this  as  the  communication  model  throughout  the  paper.)  If 
we  define 


xt  =  Vt  cos (</>c), 

(2) 

and 

Vt  =  vt  sin((/>c), 

(3) 

the  clutter  process  can  also  be  represented  as 

ct  =  Xt  cos  ( IV ct )  -  yt  sin  (uct) . 

(4) 

If  we  take  vt  to  be  a  zero  mean  process  with  a  variance  a2,  and 
assume  the  phase  noise  is  not  a  random  process,  then  this  is  the 
standard  Gaussian  process  with  the  clutter  being  independent 
of  the  time  origin.  Thus  the  Gaussian  clutter  process  is  a 
special  case  of  the  communication  model.  Similarly  the  Rician 
Clutter  model  is  just  the  Gaussian  model  with  a  scatter  ( S ) 
added  to  the  mixture.  The  log-normal  clutter  model  is  to 
take  the  communication  model  for  clutter  with  the  envelope 
probability  distribution  taken  to  Rayleigh  distributed  and  the 
phase  noise  taken  to  be  uniformly  distributed.  Other  models 
of  the  amplitude  that  have  been  used  include  the  Weibull 
representation  and  the  K-distribution.  The  problem  with  these 
approaches  to  analytical  clutter  is  they  treat  the  amplitude 
differently,  but  the  not  the  phase.  All  analytical  models  assume 
that  the  phase  noise  is  uniformly  distributed.  This  is  clearly 
not  the  case,  but  there  has  been  no  method  that  overcomes  this 
difficulty.  In  this  paper  we  propose  a  method  for  characterizing 
non-uniform  phase  noise,  so  one  can  determine  the  phase 
noise  clutter  distribution  and  then  use  this  with  a  desired 
amplitude  distribution  determine  better  analytical  models  for 
clutter.  We  note  that  the  analytical  clutter  determination  using 
the  communication  model  of  clutter  is  a  special  instance  of 
the  Rayleigh  Problem  [4],  [4],  [5]  originally  arose  in  Lord 
Rayleigh’s  investigations  into  scattering  acoustical  waves  off 
of  rough  surfaces.  The  problem  was  to  evaluate  the  PDF  of 
the  random  sums  such  as 

N  N 

X  =  ^2  Xi  =  ^2  at  cosCdi).  (5) 

i= 1  i=l 


2.  DETERMINATION  OF  PROBABILITY  DENSITY 
FUNCTIONS  USING  CHARCTERISTIC  FUNCTIONS 

The  communication  model  of  clutter  is  the  N  =  1  instance 
of  the  Rayleigh  problem.  In  many  engineering  applications, 
difficulty  occurs  when  we  have  a  change  of  variables.  We  know 
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the  distribution  of  the  original  random  variable  x,  but  we  have 
to  deal  with  a  change  of  variables  y  =  f(x),  so  we  are  trying 
to  cope  with  a  design  based  on  knowledge  of  x  rather  than 
y  which  is  what  we  actually  need.  A  method  for  finding  the 
underlying  probability  distribution  for  y  is  to  use  characteristic 
functions.  The  definition  of  a  characteristic  function  (CF)  for 
a  density  P(x)  is 

/OO 

ei9xP{x)  dx  (6) 

-OO 

which  can  also  be  interpreted  as  the  expected  value  or  average, 
(e*ex),  of  el9x  (see  [?]  for  a  thorough  discussion  of  CF’s).  The 
CF  has  the  properties: 

1)  M( 0)  =  1. 

2)  M*(-0)  =  M(6). 

3)  \M{6)\  <  1. 

4)  \M{0)\  <  M(0). 

There  are  several  methods  for  defining  the  Fourier  transform 
of  a  function  that  are  used  in  the  literature;  all  involve  the 
placement  of  the  factor  27t  that  always  occurs  when  using  the 
Fourier  transform.  We  use  the  asymmetric  form  that  is  found 
in  Bracewell,  Papoulis,  and  others.  If  we  adopt  the  Hilbert 
space  notation  for  the  inter-product  of  two  functions  (where  * 
denotes  complex  conjugate) 

/OO 

dt,  (7) 

-OO 

then  the  Fourier  transform  of  a  function  /(f)  can  be  defined 
in  Hilbert  space  notation  as 

/OO 

eiut  f(t)  dt.  (8) 

-OO 

The  inverse  Fourier  transform  of  a  function  F(jjj)  is  defined 
as 

1  1  r°° 

/(*)  =  ^  (eiut,  F(u))  =  —  J  e~iu)t  F{w)  du>.  (9) 

Also,  note  the  definition  of  the  Dirac  delta  function  in  this 
notation  is 

1  1  f°° 

S|t|=s(‘  <>») 

Note  the  definitions  and  results  we  use  are  consistent  with  the 
various  Papoulis  books  with  i  =  —j  as  the  translation  device. 

In  many  engineering  applications,  difficulty  occurs  when  we 
have  change  of  variables.  If  one  has  a  random  variable  x  with 
a  density  P(x),  one  would  often  like  to  know  the  density  of 
a  new  variable,  say  u,  that  is  a  function  of  the  old  variable 

u=  f(x).  (11) 

The  standard  method  for  solving  this  type  of  problem  is  found 
in  [7],  Assume  knowledge  of  the  characteristic  function  is 
known,  then  the  distribution  can  be  obtained  by  inversion.  The 
characteristic  function  is  the  expected  value  of  e*w  P x ■* 

Mu(u)  =  =  (e~iu^x\Px{ x))  .  (12) 


Now  if  we  substitute  this  definition  of  Mu(lo)  into  the  defin¬ 
ition  of  P(x)  gives 

pu(u)  =  ^  (eiuu,(eiuRx\Px{x)))  (13) 

which  upon  rearrangement  of  the  order  of  integration  gives 

p(u)  =  h  ((ei“Ue~iulf(x)’1) >  p*(x))  (14> 

follows  from  the  definition  of  the  delta  function  [8]  that 


d(u-f(x))  =  ^(ei“(u-px»,  l)  (15) 

P(u)  =  (S(u- f(x)),Px(x)) .  (16) 


Now  one  of  the  properties  of  the  delta  function  is  [9] 

*(»(*))  =  £«(*-*.)  5753, 


(17) 


were  the  x^s  are  the  solution  to  g{xi)  =  0,  e.g.  the  zeros  of 
g.  Now  for  our  problem,  g(x)  =  f(x)  —  u,  the  derivatives  of 
g  is  the  derivative  of  /,  so  we  have 


6(f(x)-u)  =  Yl5(x~Xi)\f'(x.)\- 

Thus  the  expression  for  the  density  function  becomes 


I  f'{Xi)\’ 


(18) 


(19) 


If  f  is  single  valued,  then  the  solution  is  particularly  simple. 

The  transformation  y  =  Rsin(ip)  a  PDF  (<, p )  is  onto 
but  not  one-to-one  over  the  interval  [— oo,  oo].  Thus  it  has  an 
infinite  number  of  zeros.  It  is  more  convenient  to  determine 
the  CF  directly,  so  the  Fourier  transformation  of  the  PDF  is 

(20) 

The  exponential  can  be  written  as 


Jn^R)ejnv>  =  eiuRsin^\  (21) 


so  the  CF  is  given  by  [10], [11] 

OO  OO 

Mv(u)  =  J2  (einV’  /M)  =  Y,  Jn(coR)F(n) 

n=— oo  n=  —  oo 

(22) 

where  where  F(n)  is  the  Fourier  transform  of  the  PDF  for  the 
angle  variable  fv(<p)  which  is  evaluated  n.  Noting  the  Bessel 
functions  can  be  rewritten  as  ( J-n(x )  =  (— )"  Jn(x)),  we  have 

OO 

M8inM  =  MuR)FM  +  Y,  Jn{ojR)S{n)  (23) 

n=  1 

where  S(n)  =  [F(n)  +  (— )nF(— n)].  If  F(n)  is  even  or  odd, 
this  formula  can  be  simplified  further.  For  the  transformation 
x  =  f?.cos(£>)  =  f?sin(£>—  f),  which  amounts  to  replacing  ip 
by  ip  —  |  in  the  exponential,  so  the  CF  is 

OO 

M™( w)  =  J0(uR)F( 0)  +  Y  Jn(uR)C(n).  (24) 

n=  1 
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where  C(n )  =  [(— )nF(ri)  +  F(— n)\.  The  expressions  for 
C(n)  and  S(n)  can  be  simplified  further  if  the  functions  are 
even  or  odd. 

The  expression  for  the  moments  of  the  characteristic  func¬ 
tion  can  be  expressed  as  the  derivative  of  the  CF, 

1  dmMP(io)\ 


(xm)  =  + 


<9u+ 


uj— 0 


Rn 


{x™)=  —  2™{S{p)/C{p)}J(irXuR) 


(25) 


(26) 


cj—0 


ae2  ^  e 


\/2<707r 
thus  F{n )  =  e_<T9"2/2 


(27) 


(,2)  = 


R 2 


1  —  e 


(28) 


<*4>  = 


d4Mv(to) 


dto 4 


^  1 

1 

II 

e~8a°  -  4e_2<T®  +  3 

00 

O 

II 

3 

For  non-zero  mean  Gaussian,  the  CF  is 


r2„2 


e-w2cr^/2e-iw6»o 


thus 


_  p-n2<Jg/2e-in80, 


F(n)  =  e~ 

Note  the  mean  is 

(x)  =  2e-<T®/2  [cos6*o]  f?, 
so  the  second  moment  is 


(29) 


(30) 


(31) 


R 


2  r 


(x2)  =  —  1  -  e  2a°  cos  (26»0) 


2  L 

while  the  fourth  moment  is  for  either  transformation 
R4  r 


(32) 


(a;4)  =  —  e  8a<>  cos(40o)  —  4e  2<Te  cos(20o)  +  6 


2)  Laplace  Distribution:  The  PDF’s  CF  is 

,,2 


_P-<AA 


<+> 


a 


thus 


=  -2^»  (35> 
ccz  +  n- 


which  is  even.  Note  the  mean  is  zero  and  the  second 
moment  is 

,,2 

11-  ' 


td2 

- 1 


a 


a * 


while  the  fourth  moment  is 


<*4>  = 


R4 


a 


-4- 


a 


To  determine  ( xn )  requires  us  know  the  n-th  derivative  of 
the  Bessel  function  since  the  characteristic  function  for  phase 
noise  always  have  Bessel  function  components.  Since  ./(> (0)  = 
1  and  J„(0)  =  0  for  n  ^  0,  the  only  terms  that  remain  after 
we  take  the  derivative  with  respect  to  u>  and  set  it  equal  to 
zero  are  those  Bessel  functions  that  have  zero  coefficients, 
e.g.  those  of  the  form  Jp_m( x)  which  are  one  when  p  =  m. 
allows  us  to  determine  the  moments  to  arbitrary  order.  Three 
examples  of  distributions  are: 

1)  Gaussian  or  Normal  Distribution:  The  PDF  of  the  zero 
mean  normal  distribution  has  a  CF  given  by 


3)  Chi  Squared  x2(rc):  The  PDF’s  CF  is 

xm/2-l 


2"*/2r(m/2) 6  X/'&{X)  ^  (1  _  2?'w)m/2  ’ 


where 


thus 


©(+)  = 

F(n)  = 

Note  that  the  mean  is 

/  \  ■  1  1 
W  =  * 


1  x  >  0 

0  a;  <  0  ’ 

1 

(1-2  in)m/2' 


1 


which  is  even.  Note  the  mean  is  zero  and  the  second 
moment  for  either  transformation  is 


(1  +  2  i)m/2  (1-2  if1'2 

and  the  second  moment  is 

1  1 


R 


(36) 


(37) 


(38) 


(39) 


(40) 


(41) 


<*2>  = 


2  - 


(1  -  4 j)m/2  (1  +  4 i)m/2 


R2.  (42) 


The  third  moment  is  zero,  and  the  fourth  moment  for 
either  transformation  is 


From  the  characteristic  function,  we  can  determine  the 
PDF  of  the  transformation  y  =  R  sin (C)  by  using  a  Fourier 
transform  identity  found  in  [12]  which  enables  us  to  determine 
the  PDF  of  the  sinusoidal  transforms  from  the  expression  for 
Mv(u>)  to  give 

jsin8  ^  _  2  [^(0)  +  En=l  ^sin(g)  {n)Tn ( |r )]  _  j/  ^ 


where 


\/(i-(*)2> 
v. <„<«)(»)  =  [F(n)  +  (-)"?(-»)]  (-<)' 


(43) 

(44) 


and  Tn(-)  are  the  Chebyshev  polynomials.  The  PDF  of  the 
transformation  y  =  Rcos(9)  is 


jcos 6(  }  _  2  [F(0)  +  J2n=l  Ucos(8){n)Tn{  Vu) 


where 


UcosW(n)  =  [(-)nF(n)  +  F(-n)\. 


(45) 

(46) 


Cl 


■  u 


(33) 

(34) 


To  determine  the  PDF  from  a  given  characteristic  function,  we 
merely  substitute  the  particular  F(n)  and  we  have  the  under¬ 
lying  PDF  for  sinusoidal  transformations  for  any  distribution. 

3.  TWO  DIMENSIONAL  PDF’s  FOR  ANALYTIC  CLUTTER 
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Since  clutter  results  from  scattering  off  a  rough  surface, 
there  are  four  possibilities  to  consider  for  either  clutter  or 
scattering: 

1)  All  the  amplitudes  are  deterministic  and  individual  ran¬ 
dom  phase  angles  are  the  same  type  of  distribution  with 
different  means  and  variances, 

2)  All  the  amplitudes  are  deterministic  and  individual  ran¬ 
dom  phase  angles  have  different  distributions, 

3)  All  the  amplitudes  are  random  with  he  same  type  of 
distribution  with  different  means  and  variances  and 
individual  random  phase  angles  are  the  same  type  of 
distribution  with  different  means  and  variances, 

4)  All  the  amplitudes  are  random  with  he  same  type  of 
distribution  with  different  means  and  variances  and 
individual  random  phase  angles  are  phase  angles  have 
different  distributions. 

Once  one  has  the  individual  z,  then  one  has  to  determine 
the  PDF  of  N  random  sums 

Z  =  zi  +  Z2  +  ■■■  +  zn-  (47) 


To  evaluate  such  sums,  it  suffices  to  calculate  them  two-fold 
at  a  time,  so  sums  such  as  2)2  =  zx  +  22,  with  distributions 
PZl(z\)  and  PZ2(z2)  (assuming  they  are  uncorrelated),  then 
the  distribution  of  2  which  is  denoted  by  Pz  (2)  is 


/OO 

PZlz2{zi,z~  21)  dzi 

-OO 

/OO 

PZl(zi)PZ2(z  —  Z\)  dz\  (uncorrelated)  (48) 

-OO 


To  compute  the  PDF  of  Z,  we  just  continue  the  process 


Z 123  —  Z12  +  +23, 


(49) 


so 

Z  =  Z123...N-1  +  zn-  (50) 

In  principle  the  problem  is  solved  provided  the  PDF  for  fz(z ) 
can  be  calculated  so  one  can  then  calculate  Pz(z).  Calculations 
of  convolution  of  Chebyshev  polynomials  are  not  difficult,  so 
we  can  calculate.all  of  these  sums. 

Note  if  one  is  solely  concerned  with  the  first  and  second 
moments  of  the  distributions,  the  individual  distributions  can 
be  used  to  compute  the  collective  means  and  standard  devi¬ 
ations  for  Z.  For  the  Gaussian  distribution,  the  mean  is  zero 
and  the  second  moment  is  (where  we  have  replaced  R  with 
the  amplitude  A) 

(x2)  =  2A2  (l  -  e~2a2^  (51) 

Thus  for  case  1  we  have 

(Z2)  =  2NA2  (l  -  e"2<T»)  ,  (52) 

for  case  2  we  have 

N 

(Z2)  =  2NA2  -2A2^e_2<7%  (53) 

for  case  3  we  have 

N 

(Z2)=2(l^e-2^)^A2,  (54) 

i— 1 


and  for  case  4  we  have 


N 

(Z2)  =  2]>>2( 


1  —  e 


—  2  <ri 


i= 1 


Note  when  the  mean  is  not  zero,  we  have 
( Z )  =  2e-0’f)/2  [cos  0o]  A 
and  the  second  moment  is 

(Z2)  =  2  [l  -  e"2^  cos  (20o) 


A2 


(55) 


(56) 


(57) 


and  we  replace  the  previous  cases  with  e  2a<>  cos  (20o)  for 
e  6  with  the  proper  index  as  necessary.  State  of  the  art  signal 
processing  requires  using  statistics  beyond  second  order  so  that 
requires  returning  to  PDF’s  of  the  scatters  and  evaluating  their 
convolutions. 

Now  we  know  the  functional  form  for  the  sinusoidal  trans¬ 
form  is  either  (43)  or  (44).  For  the  last  two  cases  we  need 
the  formula  for  the  density  function  for  the  product  of  two 
random  variables  2  =  st  which  is  given  by 

r°°  1  z  r°°  1  z 

fz(z)=  ]— 1  fst(w,-)dw=  ]— |  fs(w)ft(-)dw, 

7-00  M  W  7-00  H  w 

(58) 

when  s  and  t  are  independent.  (Note,  by  symmetry,  it  doesn’t 
matter  which  order  of  labeling,  s  or  f  one  chooses  in  the 
integral.)  It  does  matter  whether  we  use  sine  or  cosine  since 
once  one  integral  is  evaluated  the  other  can  be  evaluated 
in  terms  of  the  first,  so  the  cosine  is  considered  in  the 
remainder  of  the  paper.  The  probability  density  function  for 
the  transformation  2  =  a  cos  (j)^j  is  (fa[w)  is  the  PDF  for  a) 

UZ)  =  2  f  ^  +  ^UcosWi^Tniw)]  z_  ^ 
7T  7-1  M  V  1  -  W2  w 

(59) 

There  are  two  limiting  cases  for  fa(w )  that  are  of  particular 
interest,  a  uniformly  distributed  amplitude 

1 


/»M  = -e(i-|[[|) 

and  a  normally  distributed  amplitude 

1 


fa(w)  = 


;  exp  - 


TTCTt 


w 


2ct2 


(60) 


(61) 


For  the  uniformly  distributed  amplitude 


fz(z)  = 


m  f1  e(i- 1^1) 


7 ra 


v^ 


dw 


-1  |  to  | 

EZiUcos{e)(n)  Tn(w)Q  (1  -  |^|) 


air 


m) 


Now  the  first  integral  is  even,  and  it  can  be  evaluated  in  closed 
form  [13]  341.01 

|  b  +  \fb2  —  x2  | 
x 


1 


1 


M  y/b2  —  x 


2  = 


First  note  that  Tn(—w)  =  (—)nTn(w),  and  [14] 

In/ 2] 

Tn{x)  =  Y.  C^X^ 

k= 0 


(63) 


U.S.  Government  work  not  protected  by  U.S.  copyright 


( n )  _  2n-2k-l 


Now  the  definition  of  the  beta  function  is 


B(m,n)  =  f  xm  1(\—x)n  1dx 
Jo 


while  definition  of  the  gamma  function  is 


r  (m  +  i)  = 


B(m,n)  = 


Now  the  first  integral  in  the  PDF  can  be  evaluated  as  (k(z)  = 

(64)  2^l> 

f1  exP 

/  - - o  dw  (78) 

Jo  wy\  —  wz 

(65)  =  If-,  r(2n-l)!!n,„(n+l) 

2  ^  (2n)!!fc(z)n+i 

=  'v;(-)"^m(,))  (so) 


The  second  integral  can  be  expressed  as 


r(?n)r(n) 


r(m  +  ?i) 

The  definition  of  the  incomplete  beta  function  is 

ru 

Bu(a,b )  =  /  xa~1(l  —  x)b~1dx  Re(a)  >  0 
Jo 

while  the  complementary  incomplete  beta  function  is 

B^(a,b)  =  f  xa~1(l  —  x)b~1dx  Re(a)  >  0 

J  U 


Bu(a,b)  +  Bcu(a,b)  =  B(a,b).  (70) 

Then  the  integral  is  expressible  in  terms  of  the  complementary 
incomplete  beta  function  as 


Tn{w)  exp  (-3,^)  dw 


M  vi  —  w2 


1  (4m  —  l)!!fT(m,  k(z)) 

2“ \  (4  m)\lk(z)2m+^ 


where 


r 1  w2m—2k  —  l 


dw  =  Bca  [m  —  k  —  1, 


H(m , k(z))  =  rk(z)(l  +  n-m+ -) 

since  it  is  even.  Thus,  PDF  for  an  arbitrary  phase  noise  and 
Gaussian  amplitude  noise  follows  from  these  integrals  and  is 
given  by 

,  ,  ^  ^(0)  -,n  (2n  —  1)!!  ,  ,,  „ 

fziz)  =  — ==>(-)  R{n,k{z )) 

^V^l^o  (2n)!! 

V°°  77  to™'!  (4m-iy.i.H(m,k(z)) 

Em=i  Ucos(e){2m) 


The  PDF  for  an  arbitrary  phase  noise  and  uniform  amplitude  We  haye  shown  it  is  possible  t0  determine  the  PDF  for  specific 

noise  follows  from  these  integrals  and  is  models  of  amplitude  noise  and  arbitrary  phase  noise.  From 

.F(O)  7  _|_  yj ~2  __  a2  these  expressions,  one  can  work  out  either  the  moments  or 

fz{z)  =  - -  log  -  +  (72)the  characteristic  function  of  the  distributions.  Once  the  PDF 

’  has  been  determined  for  z  ,  sums  can  be  considered.  Let  us 

Em= 1  Ucos(e)(2m)  EL=o  c^m)Bl  (m  —  k  —  1  J)  consider  the  two  sum  case  for  each  of  these  possibilities  and 

- - - 5 - (73)  see  what  conclusions  we  can  draw.  For  the  two  sum  case  Z  = 

njr  --v. 

z\  +  Z2,  we  know  the  PDF  of  Z  is 

For  the  normally  distributed  amplitude,  the  PDF  is  roo 

,  2  x  Pz(z)  (unco^elated)  /  PZl{zi)PZa(zi  —  z)  dz\.  (83) 

2F(0)  /•lexp(-2#FjJ  j  ,  _  J~°° 

Jz\z)  —  f- — »  /  r - 5-  aw  +  Thus  our  ability  to  determine  Pz(z)  depends  on  our  ability 

7T  sJ2tto1  Jq  wyl  -  W  ,  ,  .  .  . 

v  to  evaluate  i°r  either  sine  or  cosine  integrals  of  the  form  so 

E°°  U  (6)(n)  f1  exP  y  2w2a2  )  dw  our  anzatz  for  the  general  case  applies  by  computing  the  sums 

two-fold  at  a  time.  Details  of  these  calculations  and  explicit 
TTs/Ziraf,  j-i  \w  VI  —  w  .  .  .  ,  ....  .  , 

v  determination  of  the  moments  will  be  considered  in  a  future 

The  definition  of  the  incomplete  gamma  function  is  paper. 


z  +  v  z2  —  or 


fz{z)  = 


rz(m+  1)  =  f  xme  Xdx, 
Jo 
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while  the  definition  of  the  complementary  incomplete  gamma 
function  is 


/»00 

r®(m+l)=  /  xme~xdx, 

J  Z 
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(76)  Code  B-10 

Dahlgren,  VA  22448 
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